Cancellation of Spin-Orbit Effects in Quantum Gates Based 
on the Exchange Coupling in Quantum Dots 
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We study the effect of the spin-orbit interaction on quantum gate operations based on the spin exchange 
coupling where the qubit is represented by the electron spin in a quantum dot or a similar nanostructure. 
Our main result is the exact cancellation of the spin-orbit effects in the sequence producing the quantum 
XOR gate for the ideal case where the pulse shapes of the exchange and spin-orbit interactions are iden- 
tical. For the non-ideal case, the two pulse shapes can be made almost identical and the gate error is 
strongly suppressed by two small parameters, the spin-orbit constant and the deviation of the two pulse 
shapes. We show that the dipole-dipole interaction leads only to very small errors in the XOR gate. 
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The spin 1/2 of an electron is a "natural" representa- 
tion of a quantum bit (qubit) since it comprises exactly 
two levels; there are no additional degrees of freedom into 
which the system could "leak" and thereby cause errors 
in a quantum computation. In addition to this, magneto- 
optical experiments have revealed unusually long spin co- 
herence times in doped semiconductors, exceeding 100 fis 
[jjj , thus making electron spins in semiconductors suitable 
candidates for a scalable quantum computer architecture. 
These advantages have motivated the idea of spin-based 
solid-state quantum computation using electron spins in 
coupled quantum dots pj, where the required two-spin 
coupling is provided by the Heisenberg exchange inter- 
action between the two spins in adjacent quantum dots. 
The microscopic origin of the exchange coupling lies in 
the virtual tunneling of electrons from one quantum dot 
to the other and back, and there are several external 
physical parameters (gate voltages, magnetic field, etc.) 
which can in principle be used for controlled quantum 
gate operation M. Subsequent schemes Q ||, g] 
for solid-state quantum computation rely also on the ex- 
change interaction between spins, and it has been pointed 
out that the exchange interaction alone (without single- 
spin manipulation) is in principle sufficient for universal 
quantum computation ]9|, |l^| . 

However, the two-level structure of the spin of the elec- 
tron is only approximate if one includes relativistic effects 
which lead to spin-orbit coupling fl])| . The exchange 
Hamiltonian can acquire anisotropic terms due to spin- 
orbit coupling [jl2| |l3) . For conduction band electrons in 
single GaAs dots, the spin-orbit energy is typically small 
H , however it was recently pointed out by Kavokin |TJ] 
that the spin-orbit coupling can be relevant for tunneling 
between two dots, leading to an anisotropy in the result- 
ing spin Hamiltonian, and it was suggested that it may 
lead to additional spin decoherence. Subsequently, Bon- 
esteel et al. |l5| have demonstrated that the first-order 
effect of the spin-orbit coupling during quantum gate op- 
erations can be eliminated by using time-symmetric pulse 
shapes for the coupling between the spins. 



In this paper, we present a different method for deal- 
ing with the spin-orbit interaction. Our main result is 
that the spin-orbit effects exactly cancel in the gate se- 
quence [Eq. (||)] required to produce the quantum XOR 
(CNOT) gate, provided that the pulse form for the spin- 
orbit and the exchange couplings are the same. Since 
XOR is sufficient to assemble any quantum computation 
together with single-qubit operations, this result has far- 
reaching consequences for spin-based quantum compu- 
tation with the exchange interaction; it ascertains that 
the spin-orbit coupling can be dealt with in any quan- 
tum computation. In reality, the pulse shapes for the ex- 
change and the spin-orbit coupling cannot be chosen com- 
pletely identical. Typically, however, we can to choose 
two pulse shapes which are very similar and show that our 
result still holds to a very good approximation, i.e. the ef- 
fect of the spin-orbit coupling is still strongly suppressed. 
Finally, we discuss the effect of the dipole coupling be- 
tween adjacent spins, providing another anisotropic cou- 
pling. The anisotropy due to an inhomogeneous magnetic 
field was studied in ]hJ . 

The spin-orbit coupling for a conduction-band elec- 
tron (momentum k, spin S) can be written as H so = 
h(k) • S. In two dimensions, the Rashba term (ItJ 
hi(k) = a\{k y , — k x , 0) arises from an asymmetric quan- 
tum well or from an external field. The absence of the 
inversion symmetry, e.g. in GaAs, causes a term |l6f| 
Ii2(k) = a2(—k x , k y , 0). Such a term was already shown 
to exist in |19f| . The isotropic Heisenberg coupling with 
exchange energy J and the anisotropic exchange between 
two localized spins Si and S2 (s = 1/2) are combined in 
the Hamiltonian || H(t) = J(t) (Si • S 2 +A(t)). We 
divide A(t) into asymmetric and symmetric parts 



A(t) = (3(t) ■ (Si x S a ) + 7(<)09(*) ' S i)(/3W ' S 2 ), 



(1) 
the 



where (3 = (■0i|«h(k)|V'2) is the spin-orbit field, \ip 
groundstate in dot i = 1,2, and 7 w 0(/3°). For .4 = 0, 
the quantum XOR gate can be obtained by applying 
H(t) twice, together with single-spin rotations 0, pojj , 



U a 



u. 



(2) 



2 



where U is the (unitary) time-ordered exponential U = 
Texp(— i J^ s J 2 ^ 2 H(t) dt). Here r s denotes the switching 
time, during which the spin interactions via tunneling are 
turned on. In the case .4 = 0, the Hamiltonian commutes 
with itself at different times and thus U is only a function 
of the integrated interaction strength, 



»/2 



J(t) dt , 



(3) 



-T./2 



with tp 7^ 0. In particular, we obtain the desired quantum 
gate (up to a trivial change of basis) U g — Ucpf — 
e i7rS 2 /2 U XOR er % ' nS ?l' 1 if we choose tp = n/2 (in this case, 
U is the "square-root of swap" gate ) . 

First, we study the case A ^ 0, retaining the property 
that H{t) commutes with itself at different times. This 
is the case if (3 and 7 (and thus A) are time- independent, 
i.e. if the anisotropic part of the Hamiltonian H is pro- 
portional to the isotropic exchange term. This allows us 
to fix a coordinate system in which (3 points along the z 
axis, and in which the anisotropy can be written as 



A = (i{SlS v 2 - SfS%) + SSfS. 



2 ■ 



(4) 



with S = 7/3 2 . In this basis H commutes with the z com- 
ponent S z = Sf + Sf of the total spin, [H, S z ] = 0, and 
thus I ft) an d III), being non-degenerate eigenstates of 
S z , are also eigenstates of H . Note that in their energy 
eigenvalue J(l +£)/4 there is no contribution from the 
first term in Eq. (Q) . In the S z — subspace we choose a 
basis consisting of the spin singlet \s) = (| j|) — 1 1|))/\/2 
and the triplet \t) = (|||) + |lt))/V2~ because this 
choice makes the isotropic part JSi • S2 of the Hamil- 
tonian diagonal. The complete Hamiltonian in the basis 

{|TT),k>,|*>,|U>} is 



H{t) = 



At) 



( 1+6 \ 

-10 
0-0 1 
\ 1 + 8 J 



(5) 



where we have added an irrelevant term J(l + S)/4 pro- 
portional to the unity matrix. 

Exponentiation of Eq. (||) in the S z = subspace yields 



U 



c + is/x (3s 
—(3s c — is/x 



(6) 



where c = cos(a;</?/2), s = sin(x^/2), x = y/l + /3 2 and 
where </? is defined in Eq. (||). Since exp(z7rS'|) = — icr^ in 
the S z = subspace, we find Ue lnS ^ U\s'=o = —i<J x , i-e. 
the dependence of U\s*=o on the phase y> as well as on 
the spin-orbit parameter f3 as shown in Eq. (^|) cancels 
exactly in the sequence Eq. (^|). In other words, when 
we construct the XOR gate there will be no effect of the 
time-independent anisotropic terms A in the S z = sub- 
space. By a proper choice of tp, we can also eliminate the 



effect of the anisotropy for the states | ft) and 1 1|). This 
be seen by writing down the full unitary operator 



can 



Eq. (H) using the Hamiltonian Eq. (||), 

U g = diag(ze-^ (1+5) , 1, 1, -z e -^ (1+5) ) 



(7) 



where diag(a;i, . . . , 2:4) denotes the diagonal matrix with 
diagonal entries xi, . . . ,x&. The pulse strength <p and 
the spin-orbit parameters only enter U in the S z = ±1 
subspaces. We would like U g to be the conditional phase 
flip operation Ucpf = diag(l, 1,1,-1), being equivalent 
to the XOR operation up to a basis change. Indeed, the 
condition U g = Ucpf can be fulfilled for p = 7r/2(l + 5). 

We have shown that in the case where the anisotropic 
term in the Hamiltonian is proportional to the isotropic 
term (i.e. A = const.), we can completely eliminate the 
effect of the anisotropy by a proper choice of the pulse 
strength p. In real systems, however, the anisotropic 
terms in the Hamiltonian H cannot be expected to be 
exactly proportional to J{t), i.e. A(t) is time-dependent. 
In general, both f3 and 7 depend on time. Under these 
circumstances, we cannot exactly eliminate the effect of 
the anisotropy because of the time-ordering in the defini- 
tion of U and since the Hamiltonian does, in general, not 
commute with itself at different times, [H(t), H(t')] ^ 0. 

In the following, we estimate the errors due to the 
anisotropy in the Hamiltonian in the case where -A(t) is 
only weakly time-dependent. Subsequently, we present 
a procedure which allows us to achieve exactly this sit- 
uation (i.e. a weakly time-dependent A). We write 
A(t) = Ao + AA(t), where Ao is constant, as in Eq. (Q), 
and AA(t) is the small time-dependent deviation from 
Ao- The Hamiltonian is written as the sum H(t) = 
H Q (t) + H'(t) where H (t) is given by Eq. (g) and 

H'{t) = J(t)AA(t) (8) 
= J(t) (Af3(t) ■ (Si x S 2 ) + A 7 (i)(/3 • Si)(/3 • S 2 )) . 

Note that in the symmetric part we have already omit- 
ted terms which are of order A/3A7. This Hamilto- 
nian generates a unitary time evolution U = Uq + AU, 
where Uo — Texp(— i fj'J^Hoitydt) is the contribu- 
tion due to Hq. The explicit form of AU is rather 
complicated; however, we are only interested in esti- 
mating the gate error AU g = U g — Ucpf caused by 
H'(t) (note that AU g is not unitary). For this pur- 
pose we work in the interaction picture with respect 
to H (t), where Ui = U^ U = Texp(-i J^ 2 /2 H'j(t) dt) 

and Hj = UqH'Uo. In this representation, the devi- 
ation AU from the "ideal" time evolution U becomes 
AUi =Ui-l = -i f_ 3 j 2 2 H'j(t) dt + 0{H'j 2 ). The norm 

of the gate error ||AJ7 S || = max^^i \J (ip\ AUg"AU g \ip) 
(to lowest order in H') can now be estimated as follows, 

\\AU g \\ < 2\\AU\\=2\\U AU I \\=2\\AU I \\ 

< 2t s max \\H'(t)\\ = 2A , (9) 

\t\<T s /2 
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where the first equality comes from Eq. (|2j) and the uni- 
tarity of the involved quantum gates. Using /? < 1, we 
approximate ||iJ'(i)|| < \J(t)A(3(t)\/2, since the second 
term in H' is 0(/3 2 ). We use Af3(t) = (3{t) - /3 to write 

A = ^ max ^(^-l) , (10) 

2 \t\<r,/2 J V A) / 

where Jo denotes the average exchange coupling, Jo = 
f/r s 7^ 0. Note that the position of the dots is fixed 
during the switching process, thus (3(t)/\/3(t)\ — const. 
0. For the XOR gate, tp ss tt/2. The error probability 
for the described gate operation can now be estimated as 
e = ||A* out || 2 < \\AU g ^ in \\ 2 < \\AU g \\ 2 < 4A 2 . 

In order to obtain an estimate for A, we consider the 
case of coupled quantum dots in a 2DEG. For parabolic 
confinement potential V(r) — mujr 2 /2, the ground-state 
orbitals are ip(r) — (ira 2 ^ 1 exp(— r 2 /2a 2 B ), where a B — 
^Jfi/muj is the effective Bohr radius of the electronic or- 
bitals and m is the effective electron mass. If two such 
quantum dots containing one electron each are separated 
by a distance 2a, the exchange coupling between the spins 
of the electrons at zero magnetic field is given by || 

'<*«> " Sife - 1 » + T< 1 + ^) ■ 

(ii) 

where Iq is the zeroth order Bessel function, d= a/a% the 
dimensionless ratio between the half-distance a and the 
effective Bohr radius a° B = y/h/mujo, c characterizes the 
strength of the (bare) Coulomb interaction (Tlujo — 6 meV 
and c = 1.71 in our numerical example), and q — lo/loq 
is the strength of the confinement uj in units of its min- 
imum value ujq- Following E3, we find for both hi and 
h 2 that b(d,q) = \J(d,q)f3(d,q)\ — b ^/qdexp(—2qd 2 ), 
where bo — a,i/a g , i — 1,2. For Ii2 in a 5 nm wide [100] 
GaAs quantum well a 2 ~ 2 meVnm, or (3 w 0.02 at 
d = q = 1. In Fig. [IJ we plot J(d, q) and b(d, q). 

The switching process can be modeled e.g. by a time- 
dependent distance d between the dots or by a time- 
dependent confinement strength q. Here, we choose the 
latter possibility and use a pulse q(t) = uj(t)/ujQ — 
cosh 2 (at/r s ), where we choose a = 4. This pulse shape 
is suited for adiabatic switching ||, |l| and leads to a 
pulsed exchange interaction J(t) = J(d,q(t)) and spin- 
orbit field b(t) = b(d,q(t)), where -r s /2 < t < t s /2. The 
pulse shapes of the resulting exchange coupling J(t) and 
spin-orbit field b(t) are plotted in Fig. ^. In our exam- 
ple, the switching time amounts to r s = tt/2Jq sa 140 ps. 
Note that a pulsed switching by electrostatic lowering of 
the tunneling barrier between the dots or by applying a 
magnetic field results in a very similar time-dependence 
of J and 6 and a similar analysis could also be done in 
these cases. If the pulse shapes of b and J were identi- 
cal, then the effect of the spin-orbit coupling in the XOR 
gate could be eliminated exactly (as explained above). 
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FIG. 1: The exchange coupling J (dashed line) in units of 
hujo and the spin-orbit field b = /3J (solid line) in units of 
bo for two electron spins located in adjacent quantum dots as 
a function of the dimensionless parameter q — uj/uiq at fixed 
interdot distance d — 1 (Inset: as a function of d at fixed 
q = 1), where Tioj is the single-dot confinement energy, ojo 
is the (fixed) minimum value of u, and bo is the spin-orbit 
parameter. For this plot, Tiwo = 6meV and c = I.7L 

The optimal choice of f3o (i.e. the one which minimizes 
A) in our numerical example turns out to be (3q ~ (3(t = 
0.1t s ) and from Eq. @ we find A « 7 - f0~ 3 . There- 
fore, the gate errors occur at a rate e < 4A 2 w 2 • f 0~ 4 
which is around the currently known threshold for fault 
tolerant quantum computation |22[ and could therefore 
be corrected by quantum error correction. Note that in 
cases where the error e is too large for quantum error cor- 
rection, it can be further reduced at the cost of a slower 
gate operation. This can be achieved by designing pulses 
with smaller intensity, where there is a long period of 
constant A between the rise and fall of the pulse. 

We finally include the dipole-dipole interaction 
Hd = 77 (Si • S2 — 3(Si • a)), being another source for 
anisotropic coupling among the spins Si and S2, into our 
discussion. Here a denotes the unit vector pointing from 
the center of one to that of the other dot. The coupling 
parameter 77 = ^g 2 fig / An (2a) 3 is typically much smaller 
than the spin-orbit energy, for g = 2 and a = 20 nm we 
obtain 77 sa 3 • 10~ 12 eV, corresponding to a dipole field of 
Bd = v/^B ~ 0.5 mG. Nevertheless, we show here that 
in cases where the dipole interaction matters (e.g. if g is 
very large) , it can again be dealt with by using the meth- 
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FIG. 2: Pulse form of the exchange coupling J(t) (dashed 
line) and the spin-orbit field b(t) = J(t)(3(t) (solid line) for a 
simple model involving two coupled quantum dots which are 
coupled and decoupled with the time-dependent confinement 
strength q(t) = u(t)/ujo. For the distance between the dots 
we choose d = 1. The choice of the vertical scaling of the two 
pulses in this graph is such that the deviations of one pulse 
from the other are (approximately) minimal. 
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ods described above. An essential difference between the 
dipole and the spin-orbit interactions is that the dipole 
interaction between two spins located in adjacent quan- 
tum dots with fixed distance cannot be changed by apply- 
ing gate voltages or static magnetic fields and therefore 
remains constant during the entire switching process and 
in the "idle" time between the switching. In the follow- 
ing, we assume for simplicity that A is independent of 



time, as in the first part of our discussion. In addition to 
this, we assume for the moment that J is also constant. 

When investigating the combined effect of the spin- 
orbit and dipole couplings, we will restrict ourselves to 
the two cases hi and h 2 . For h 2 , the spin-orbit field 
is parallel to the interdot coupling direction (3 || a and 
thus the second term in Hd has the same form as the 
symmetric term in Eq. (0). We set £ = 77/ J and find 



J 



( 1 + 5 - 2£ 



V 





-0 







i/3 

i + e 








1 + 5 - 2^ j 



(12) 



Using Eqs. (£|) and ([|) with 92 = 7r/2(l + 5 — 2£), we exactly obtain U g = Ucpf, i-e., the combined effect of the 
spin-orbit and dipole coupling is eliminated. For hi, we find /31a. Choosing a along the rc-axis we obtain 

/ 1 + -3£/2 ^ 



H= J - 
2 






V -3£/2 



-0 







i(3 





(13) 



i+SHJ 



r 



Setting <p = 7r/2(l + 5 + £), we obtain again U g = Ucpf, 
therefore it is possible to eliminate the spin-orbit and 
dipole coupling effects also in this case. 

In principle, the analysis for time-dependent exchange 
and spin-orbit coupling can be repeated including the 
dipole interaction. However, the dipole interaction can- 
not easily be switched on and off, and therefore, the 
"pulse shape" of the dipole interaction is a constant, i.e. 
very different from those of the exchange and spin-orbit 
couplings. Nevertheless, since the dipole interaction is 
usually very small, we can still use Eq. (||) to obtain a 
reasonable upper bound on the error by setting H ' = Hd- 
We obtain = t s t\ = \ip\rj/Jo which for typical num- 
bers (as above, Jo « meV) is tiny, w 10 -9 . The 
error td — 4A^ = 4(t s t]) 2 caused by the dipole interac- 
tion is therefore negligible in typical situations, and we 
only have to take it into account if for some reason (e.g. 
large g) the dipole interaction becomes unusually large. 

We conclude that while the spin-orbit interaction can 
cause weak decoherence in the combination with phonons 
p3| , its direct effect on quantum gate operations nearly 
cancels if the pulse shapes of the exchange and spin-orbit 
couplings are as similar as possible. We have shown that 
in a typical case involving two tunnel-coupled quantum 
dots this is easily achievable. A simple estimate shows 
that the dipole interaction between the spins is usually 
much smaller than the spin-orbit interaction and can be 
neglected. Nevertheless, we have shown that in cases 
where the dipole effects are unusually large, the combined 
effect of spin-orbit and dipole coupling can be corrected. 

We acknowledge support from the Swiss NSF and from 
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